Lowering the critical temperature with eight-quark interactions 
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It is shown that eight-quark interactions, which are needed to stabilize the ground state of the 
combined three flavor Nambu - Jona-Lasinio and 't Hooft Lagrangians, play also an important role 
in determining the critical temperature at which transitions occur from the dynamically broken 
chiral phase to the symmetric phase. 
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For a long time already, a phase transition which al- 
most certainly occurs in hadronic matter at finite tem- 
perature is at the centre of attention of many researches 
and reviews [H 0, [H, HI ■ It may be very important in the 
context of heavy-ion collisions and the evolution of the 
early Universe. At present, experimental facilities with 
ultra-relativistic heavy-ion reactions aiming at studying 
the QCD phase diagram have triggered even more intense 
studies from the theoretical side, ranging from effective 
models to calculations on the lattice. Recent lattice re- 
sults predict for vanishing chemical potential and mas- 
sive quarks a non singular crossover for the QCD tran- 
sition which lead to different critical temperatures, de- 
pending on the considered observables. The transition 
related to the renormalized chiral susceptibility yields 
T c = 151 (3) (3) MeV, whereas for instance the inclusion 
of the Polyakov loop increases this value by ~ 25 MeV 

An effective chiral model which has been widely used to 
extract the low energy characteristics of hadrons involv- 
ing the u, d, s quarks is the well known Nambu - Jona- 
Lasinio model p, generalized to the SU(3)l x SU(3)r 
chiral symmetry and to include the U(1)a symmetry 
breaking of low energy QCD in form of the 2Nf determi- 
nantal interaction of 't Hooft (Nf denoting the number 
of flavors) 0- This combined Lagrangian (we denote it 
by NJLH) has been first analyzed in Q and @ to study 
dynamical breakdown of chiral symmetry and related me- 
son spectra in the vacuum, and since then in further nu- 
merous applications [l(| EH EH EH, E3 • I n more recent 
calculations the model is used to obtain the thermody- 
namic properties of hadrons and restoration of chiral and 
U(1)a symmetries, both at zero and finite chemical po- 
tential 15], and extended to include the colored diquark 
channels, of importance at large baryonic densities [lt| . 
Calculations performed with the NJLH at zero chemical 
potential yield a critical temperature for the light quarks 
around T c ~ 200 MeV, where the condensate undergoes 
rapid changes. This value clearly overestimates the lat- 
tice prediction. 

In the present letter, we argue that by extending the 
NJLH model to include interactions of eight quarks, the 
critical temperature is lowered considerably. The eight 



quark vertices have been introduced first in [171 ] to stabi- 
lize the scalar effective potential derived in the presence 
of (2Nf = 6) 't Hooft interactions [l8| and have been 
discussed at length in [l9j . 

We consider the following Lagrangian density 



£eff = 9(*7 P ^ - m)q + C iq + C 6q + C S q + ■ ■ 



(1) 



Quarks q have color (N c — 3) and flavor (V/ = 3) indices 
which are suppressed here. We suppose that multi-quark 
interactions C^ q , C§ q , Cs q are 



G 



C 6q = K(det qP L q + det qP R q), 
r - j- 

where the flavor space matrices A a , a — 0,1, ... , 

normalized such that tr(A a A/,) = 2<5 a &, and Ao = 

A Q at a ^ are the standard SU(3) Gell-Mann matrices. 
Chiral projectors are Pl,r — (1 T75)/2, the determinant 
in the 't Hooft Lagrangian is over flavor indices. The 
eight-quark interactions are given by 



(2) 
(3) 
(4) 

8, are 
1, 



(l) 2 
Hiq = 8 ffl MiPRlmMmPLqi)] 



(5) 



C 



(2) 
8„ 



16g2(qiPRq m )(q m PLqj)(qjPRqk)(qkPLqi)- (6) 



Here the summation over repeated flavor indices is as- 
sumed. One sees that the most general spin zero eight- 
quark interaction is composed of two chiral invariant 
combinations of products of four quark bilinears with 
couplings 3i,<?2, of which one, C% q \ is OZI violating. 

The model has eight parameters: the couplings of 
multi-quark interactions G, k, gi, (72; three masses of cur- 
rent quarks, m.j, and the cutoff A (the model being not 
renormalizable) . 

Since there is some freedom in the model parameter 
choices that in the realistic case with SU(2)i x U(1)y 
flavor symmetry and non zero current quark masses, 
m u = rrid 7^ m s , lead to qualitatively similar spectra 
for the low lying pseudoscalars and scalars, we chose as 
input several sets of parameters at T = 0. 
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We put furthermore all current quark masses to zero 
to have clean signatures of the phase transitions. This 
will not affect the trends we are about to discuss in con- 
nection with the pivot temperature values Ti (see below) 
obtained with and without the stabilizing eight-quark in- 
teractions. On the other hand, this simplifies essentially 
our analysis. 

Let us recall that at T = the effective potential of 
the model as a function of the constituent quark mass, 
M, in the SU(3) limit is [H 



V(M) = 



16 

3iV e 

16tt 2 



12G + Kh 



27 



ph' 



MVo(AP) +A 4 ln 1 + 



M 2 
A 2 



(7) 



with A being an ultraviolet covariant cutoff in the quark 
one-loop diagrams, p = g\ + ^g2i and 



J (M 2 ) 



9i 

A 2 - M 2 In I 1 



A2_ 

M 2 



(8) 



The function h(M) is a real solution of the stationary 
phase equation related to the integration over bosonic 
auxiliary variables (for details see e.g. [l7l ]) 



M ■ 



Gh + ^-h 2 
16 



\ph* = 0. 



(9) 



This equation yields a one-to-one real-valued mapping 
M — > h in the parameter region 



p > 0, G > 



1 



24 



(10) 



which defines the stability conditions in the SU(3) limit 
of the effective potential. 

The gap equation, obtained from the extremum con- 
dition dV/dM = 0, relates further the function h(M) to 
the quark condensate as 



h(M) 



1k. 

2ir 2 



MJ (M 2 ) = 2(0|gg|0>. (11) 



The curvature of the effective potential V(M) at the 
origin, M = 0, in the chiral limit is determined at T = 
by the value of 



N, 
2n 



2 GA > 



(12) 



with t > 1, r = 1 and r < 1 indicating a maximum, a 
saddle point and a minimum, respectively. 

At finite temperature, after introducing the standard 
non-covariant three-momentum cutoff A for the quark 
loops, see e.g. [2(|, one obtains for the effective potential 



V T (M) 



It 
16 



12G + nh 



27 



ph' 



(13) 



/ <V 



2Tln 1 + exp 
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FIG. 1: The function F(t) related with the curvature of the 
effective potential at the origin, in dimensionless units t = 
T/X. 



with E(p) = \Jp 2 + M 2 and p denoting the magnitude 
of the 3-momentum. 

The curvature at the origin is now a temperature de- 
pendent function 



27V T 
r(T ) = ^GX 2 F(t), t=£, 

7T^ A 



1 t 2 



(14) 



F W = -i + 7T^ + 2<ln C 1 + e ') + 2 * 2i 2[-e*], (15) 

2 67T Z 



where 



fc=i 



(16) 



denotes the polylogarithm function. 

The function F(t) is depicted in fig. 1. To guide the 
discussion we show in fig. 2 the temperature dependence 
of the effective potential for the parameter set (e) of Table 
1, for T = (dashed line) and three "critical" values 
TaiTd,T s , introduced as follows. 

In the zero temperature limit the curvature reduces 
to t(0) = N C GX 2 /n 2 . From this we see that if we 
wish to ajust our zero of temperature curvature so that 
t(0) = r we must choose A 2 = 2A 2 . The function F (t) 
decreases monotonically to zero as T — > oo, starting from 
F(0) = 1/2 and therefore, if the curvature t > 1 at 
T = 0, the system will undergo at some critical tempera- 
ture T — T a a change of curvature at the origin, starting 
from a saddlepoint at T a and entering a region of local 
stability near this point. Increasing further the tempera- 
ture the system undergoes changes reminiscent of a first 
order transition, passing through a configuration with 
two degenerate minima at T = until at some higher 
temperature T s a further saddlepoint appears, and for 
T > T s only the symmetric phase survives. Obviously if 
t < 1 the system is at T = already in the coexistence 
or even in the symmetric phase. 

In the following we select a few parameter sets that 
have been obtained previously with and without inclu- 
sion of the eight-quark interactions to determine the mass 
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FIG. 2: The effective potential V (in units (10 GeV)" 4 ) for 
the parameter set (e) as function of M (in GeV). Full curves: 
lowest for T a , middle for Td and upper for T s . Dashed line: 
T=0 case, the physical minimum occurs outside the plotting 
range. 

spectra of low lying pseudoscalar and scalar nonets, re- 
lated weak decay constants and quark condensates. In 
Table 1, the model parameters A, G, k, gx, g2, are taken 
from [2lJ, set (a), and from (l9j ]. sets (d,e,f), at T = 0. 
This input is kept fixed in the finite temperature evolu- 
tion of the effective potential. One sees that the eight- 
quark interactions present in sets (e,f) reduce consid- 
erably the temperatures T ai Td,T s , as compared to the 
cases with g\ = g% = 0. In the last set (f) the curva- 
ture t is slightly below 1, so the system has already a 
minimum at the origin at T = 0. 



TABLE I. Parameters of the model at T=0: G (GeV -2 ), A 
(MeV), k (GeV -5 ), gi, g 2 (GeV -8 ). Indicated are also the 
temperatures T a ,Td,T 3 (MeV)(see text). 



Sets 


A 


G 


— K 


9i 


92 


T a 


T d 


T s 


a 


820 


13.5 


1300 








192 


202 


204 


d 


839 


12.16 


1082 








174 


183 


185 


e 


839 


11.28 


1083 


1500 


327 


143 


161 


163 


f 


839 


8.92 


1083 


6000 


327 




111 


135 



Let us discuss now the reasons for the obtained de- 
crease of T c . The underlying mechanism is quite simple. 
The critical temperature is correlated with the value of 
t(0): the closer r(0) is to 1, the lower T c . The eight- 
quark interactions decrease the effective value of GA 2 (to 
see this one should consider the mass spectrum of meson 
states [Hj]), and therefore lower the transition tempera- 
ture. In the parameter sets shown we have: r(0) = 1.4 
(a), r(0) = 1.3 (d), r(0) = 1.2 (e) and r(0) = 0.95 (f). 

The lowering of the critical temperature T c has also 
been recently observed in the Nj = 2 case with eight- 
quark interactions 22]. 

We conclude that the presence of the eight-quark inter- 
actions in the hadronic vacuum results in an observable 
effect which deserves a more detailed study. The argu- 
ments presented above lead to the expectation that the 
tendency will be the same when one considers the real- 
istic case with physical current quark masses. Work in 
this direction is in progress. 
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